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New contributions to gluon poles in direct photon production 
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We consider the direct photon production in two hadron collision with one of the hadrons being 
transversely polarized. By using the contour gauge for gluon fields, we hnd that there are new 
twlst-3 terms present In the hadron tensor of the considering process in addition to the standard 
twist-3 terms. In this work, we demonstrate that the significance of these new terms are two-fold: 
first, they are crucial to both the QED and QCD gauge invariance and, second, their contributions 
to the hadron tensor are at least the same as those from the standard ones. We also study the 
resuting effects which are responsible for the universality breaking of the corresponding twist-3 
parton distributions. 

PACS numbers: 13.40.-f,12.38.Bx,12.38.Lg 


I. INTRODUCTION 

The problem of the electromagnetic (QED) gauge invariance in the deeply virtual Compton scattering (DVCS) 
and similar exclusive processes has intensively been discussed during the last few years (for example, see [T]- 
[^). This development explored the similarity with the earlier studied inclusive spin-dependent processes [5]. 
The gauge invariance of hard process amplitudes is ensured by twist three contributions and by the use of the 
equations of motion that provide a possibility to exclude the three-particle (quark-gluon) correlators from the 
amplitude. So that, after combining with the two-particle correlator contributions, one gets the gauge invariant 
expressions for the physical amplitudes or, in the case of lepton-hadron processes, for the corresponding hadron 
tensors |S] . This scheme was originally developed in the case of the particular inclusive processes with transverse 
polarized hadrons, like structure function g 2 in DIS [6] and Single Spin Asymmetry (SSA) [7] due to the soft 
quark (fermionic poles [S]). Also, the QCD gauge invariance of the so-called gluonic poles contributions [3] has 
been a subject of studies in uni where the methods that are used rely on the Wilson exponentials [IIHIl- 

We have shown in our recent work m that, in order to ensure the QED gauge invariance of the transverse 
polarized Drell-Yan (DY) hadron tensor, it is mandatory to include the contribution from the extra diagram 
originating from the non-trivial imaginary part of the corresponding twist-3 function {xi,X 2 )- Before this 
study m, the function {xi,X 2 ) has been argued as real, and the imaginary part of amplitude was ensured 
by means of a specially introduced “propagator” ^ in the hard part of the hadron tensor m- However, we have 
explained in |15| that the i?''^-function does, in fact, have an imaginary part, and the existence of this imaginary 
part can be realized with a help of the contour gauge. Moreover, the H^'^-function with the complex prescription 
induces a new contribution to the hadron tensor. As has been stressed, this extra contribution leads to the 
amplification of corresponding tensor by a factor of 2. This finding of ours has independently been confirmed 
in m by using of a different approach. Besides, from the point of view of phenomenology, the corresponding 
SSAs in the DY process and the role of gluon pole contributions have previously been discussed in . 

In the present paper, we extend our approach used in m to the case of the direct photon production in 
two hadron collision where one hadron is transversely polarized. We derive the hadron tensor for this process 
and study the effects which lead to the soft breaking of factorization (or the universality breaking) through 
the QED and QCD gauge invariance. In a similar manner as in m, the special role is played by the contour 
gauge for gluon fields. We demonstrate that the prescriptions for the gluonic poles in the twist 3 correlators are 
dictated by the prescriptions in the corresponding hard parts. We argue that the prescriptions in the gluonic 
pole contributions differ from each other depending upon the initial or final state interactions of the diagrams 
under consideration. Moreover, the different prescriptions are needed to ensure the QCD gauge invariance. We 
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treat this situation as a breaking of the universality condition resulting in factorization soft breaking. The extra 
diagram contributions, which naively do not have an imaginary phase, is discussed in detail. 

In the paper, we also show that the new (“non-standard") terms do contribute to the hadron tensor exactly 
as the “standard" terms known previously. This is exactly similar to the case of Drell-Yan process studied in 

na¬ 


il. GETTING STARTED: CASE OF DRELL-YAN PROCESS 

For pedagogical reasons, we remind briefly our findings for the Drell-Yan (DY) process where one of hadrons 
possesses the transverse polarization, see [15] for all details. As usual, the Drell-Yan process with the transversely 
polarized nucleon is defined as 

+ n{p 2 ) ^ i*{q) + x{Px) ^ eih) + iih) + x{Px) , (1) 

where the virtual photon producing the lepton pair (/i + I 2 = (?) has a large mass squared (<?^ = Q^) while 
the transverse momenta are small and integrated out. The dominant light-cone directions for the DY process 
(Fig. Ir are defined as pi « n* Q/{xb'/2) and p 2 « nQ/{ ubV^) with the dimensionless light-cone vectors 
n* ={n*~^, 0“, Oj^) and = (O’*', n“, Oj^). 

Since we deal with a large in the process under consideration, it is possible to apply the factorization 
theorem to get the corresponding hadron tensor factorized in the form of convolution: 

Hadron tensor = {Hard part (pQCD)} 0 (Soft part (npQCD)} . (2) 

Usually both the hard and soft parts in Eqn. 0 are independent of each other, UV- and IR-renormalizable. 
Moreover, various parton distributions which parametrize the soft part have to manifest the universality prop¬ 
erty. 

Based on DY-process, it is convenient to study the role of twist 3 by exploring of different kinds of asymmetries, 
for instance, the left-right asymmetry. This left-right asymmetry means the transverse momenta of the leptons 
are correlated with the direction S x where implies the transverse polarization vector of the nucleon and 
is a beam direction [S^- Generally speaking, any single spin asymmetries (SSAs) can be presented in the 
following symbolical form (at this moment, the exact expression for SSA is irrelevant): 

A , (3) 

where is an unpolarized leptonic tensor and, consequently, has only the real part; Pd^v stands for the 
hadronic tensor which is also real. Since one of hadrons is transversely polarized, the corresponding matrix 
element which forms the soft part of hadron tensor reads m 

(pi,5^|^(Ain)7+5^?(A2n)V'(0)|5^,Pi) =*e“+^"-(piP2) J dxidxa ^ 2 ), (4) 

where the light-cone vector fi is a dimensionful analog of vector n. Therefore, in order to provide for the 
condition, G ?fte, the complex i in the r.h.s. of (|^ has to be compensated either (a) by the complexness of 
the hard part (this is the standard contribution): 

~ ?9m[Hard] ® l^{p,, A)\St,Pi) ~ , (5) 

or (b) by the complexness of the soft part 

- Hard® |(pi,S'T|C>('0,V’,^)l>S'r,Pi) ^ ? 3m [B^]| . (6) 

In Eqns. ^ and ~ and -0, A) are the shorthands for the Fourier transformation and the corresponding 
quark-gluon operator, respectively (Eqn. @)- In general, the hard parts in Eqns. (|^ and ([^ differ from each 
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a) b) 

FIG. 1: The Feynman diagrams which contribute to the polarized Drell-Yan hadron tensor. 


other. For instance, the hard part of the diagram presented in Fig. [^a) contains the quark propagator in 
contrast to the hard part of the diagram in Fig. [l[b) 

However, in the previous studies (for example. ^lOl llTl l26l [27] h a; 2 )-function has been assumed to be 

a purely real function: 

B^{xi,X2) = — ^ —T(xi,X2) (7) 

Xi — X2 


with the function , 0 : 2 ) € 3?e which parametrizes the corresponding projection of {ip Gap ip)- Therefore, the 
scenario (b) (Eqn. ([^) will never be realized if the a; 2 )-function has the form as in Eqn. Q. As a 

result, the QED gauge invariance of the DY hadron tensor is in question. Indeed, having analyzed the hard 
subprocess in the context of the QED gauge invariance, we can conclude that only the sum of two diagrams in 
Fig. n (a) and (b) ensures the QED gauge invariance of the hadron tensor. As shown in [T5], the contribution 
of is associated with the Feynman diagram in Figj^a) while the contribution of is generated by the 
diagram presented in Fig. Bb)- 

Thus, we may infer that the real function {xi, X 2 ) in the form as in Eqn. (j^ leads to the violation of the 
QED gauge invariance of the hadron tensor. 

To solve this discrepancy, it is instructive to remind the reason which leads to the representation in Eqn. 
Q. The conclusion that is a real function has come from the ambiguity in the solutions of the differential 
equation (provided that A+ = 0): 


d+ = G+“. 


Indeed, the formal solutions of Eqn. 


A^iz) = 


I have the following forms: 

00 

J duj~9{z~ — + A^{—oo) 

— 00 

00 

— J du}~9{u}~ — z~)G^^{u}~) + A^{oo). 


( 8 ) 


( 9 ) 

( 10 ) 


From the first glance, Eqns. ([^ and (101 seem to be equivalent each other. However, as we will see below, this 
is not true. 

Now, if we insert Eqns. and ( (To| into Eqn. and use the following parametrization 

{pi,S^\iP{Xifi)-f+n,Gip^{X2h)iP{0)\S^,Pi) =e^+^"-{P 1 P 2 ) ^ dxidx2 T(xi, X 2 ), (11) 
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The corresponding 5-functions appeared in the hadron tensor and expressed the momentum conservation low should be also 
referred to the hard parts. This statement was argued in m in the context of the so-called factorization links. 
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we get the following representations: 


B^{Xi,X2) - 5(xi X2)PX(-oo)(a;i) + ^ J. , 

^ ' xi — X2 + le 

(12) 


(13) 


respectively. In Eqns. (12) and (13), the corresponding prescriptions ±ie arise from the integral representation 
for the theta-function: 


e{±x) 


±i 

2ti 



^—ikx 

k zLie 


— OO 


(14) 


Further, if we suppose that the representations (121 and (13) are equivalent to each other, we can calculate 
the plus and minus combinations of (121 and (13) resulting in 


B^{xi,X2) 


^B^{xi,X 2) + ^B^ixi,X2) 

^S{xi - a:2){5l(_oo)(a;i) + B^(+oo)ixi)] + ^ 


(15) 


and 


0 = B^ {xi,X2) — B^ {xi,X2) (16) 

= 5(a;i - a;2)|BX(+^)(a;i) - SX(_^)(xi)| - 2iTTS{xi - X2)T{xi,X2) ■ 


The ambiguity in the solutions of ©((EH) and ( [T^ ) ultimately gives us the standard representation ([^ with the 
real function B'^ provided the asymmetric boundary condition for gluons is given by B^^^~^{x) = ~B'^^_^-^{x). 

In fact, the representations (12) and ( [I^ are not equivalent ones m- To see that it is necessary to remember 
that all axial-type gauges are the particular cases of the most general contour gauge (see Appendix and the 
text below for details). Using the contour gauge conception, one can easily check that the representation (12) 
belongs to the gauge [a:, —oo] = 1, while the representation (13) belongs to the gauge [-|- (X), x] = 1. T herefore, 
one has no any reason to believe that (12) and (13) are equivalent ones (for details, see (|Ag, (|A:i^),he. 


{Eq.([^ ^B^(a;i,X2)} ^ {Br(xi,X2) ^ Eq.([^} . 


(17) 


Roughly speaking, it resembles the trivial situation where two different vectors have the same projection on the 
certain direction. In this context, the well-known axil gauge A+ = 0 can be treated as some type of “projection" 
which corresponds to two different ‘Vectors" represented by two different contour gauges (see Appendix [A|). 

We now consider the QED gauge invariant hadron tensor for the DY-process with the transverse polarization. 
In order to understand which contour gauge we need to deal with, before imposing the condition A+ = 0, it 
is necessary to take into account the contributions of (pi, '01‘S'^jPi) in the standard hadron tensor 

(see Fig. [^a)). Based on the analysis of the y-structure of this diagram as shown in [IS], the Feynman causal 
prescription in the quark propagator (see, eqns. (2), (6) and (8) of [TS]) uniquely leads to the Wilson line in 
the quark-gluon correlator: 


U 

[—c»“, 0“] = Pexp| — ip J dz~ {0, z~, 6 t)'^ 


(18) 


Eqn. (18) suggests that we have to use the contour gauge [—oo“, 0“] = 1. Therefore, the contour gauge defined 
by [— (X)“, 0“] = 1 destroys the ambiguity, and the function B^{xi, X 2 ) has to be described by the following 
representation (see, (A.12) and m for details) 


B+(xi,X2) 


T{xi,X 2) 
xi — X2 + ie 


(19) 
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As a result, the diagram represented in Fig. [^b) (or the hadron tensor §) does contribute to the hadron 
tensor, and this together with the first diagram represented in Figj^a) form the gauge invariant (GI) hadron 
tensor (see, Eqn. (34) of fTH]!: 


Wgi = Ww 


+ = 


( 2 ) 


=- 2 




PlP2 


[xBPi 


2 /BP 2 ) 


q{yB)T{XB,XB) 


( 20 ) 


From Fig. we can also realize that the representation of {xi, X 2 ) with the complex prescription +* e in the 
gluonic pole, see (A.121, corresponds to the initial state interaction (ISI) with respect to the hard subprocess. 
We want to stress that in the case of DY-process we deal with the initial state interaction only as opposed to 
the case of the direct photon production which is studied below. 

For the DY-process, the initial state interaction generates —2£+A:^ +ie (see, the diagram in Fig.j^a)) in the 
quark propagator which, in turn, leads to (i) the contour gauge [—oo“, 0“] = 1 and, then, to (ii) the function 
with the certain complex prescription (A.12). The latter ensures the QED gauge invariance for the hadron 
tensor. Schematically, the mentioned logical chain can be presented as 


ISI 


-£+ + ie 


gauge [z , -00 


= 1 ^ 


T{xi,X2) 


xi — X2 + le 


GI. 


( 21 ) 


We can see that the prescription in the quark propagator of the hard part gives information on the contour 
gauge for gluons from the soft part. In other words the hard and soft parts are not fully independent of each 
other. Despite of this, the DY hadron tensor has formally been factorized with the mathematical convolution, 
and the parton distributions, such as the twist-3 function B^{xi,X 2 ), still satisfy the universality condition. In 
contrast to the DY-process, as we will see in the next section, the direct photon production tensor is built with 
the functions B^{xi,X 2 ) that will not manifest the universality. 

We will refer to a soft breaking of factorization when the factorization procedure results in the mathematical 
convolution between the finite hard and soft parts, but there is no universality for the soft functions or the hard 
and soft parts are not totally independent. 


To conclude this section, let us note that the situation with the gauge invariance for the DY-process (211 


is very similar to what has been discussed for the vector meson electroproduction in m where the final state 
interaction also predetermined the correct prescription for the spurious gluon pole for the validity of QCD 
factorization. 


III. HADRON TENSOR OF THE DIRECT PHOTON PRODUCTION I: KINEMATICS AND 

GAUGE INVARIANCE 

A. Kinematics 


In this section, we study the two hadron collisions, where one of the hadrons possesses the transverse polar¬ 
ization, which produce the direct photon in the final state in: 

(Pi) + V(P 2 ) ^ 7(9) + q{k) + X{Px) . (22) 

The gluonic poles are being manifested in this process in the similar way as for the Drell-Yan process [55]. We 
perform our calculations within a coZZznear factorization, and it is convenient (see, e.g., |34| 1 to fix the dominant 
light-cone directions as 

Pi = n* , p 2 = n , with 

n; = (1/v^, Or, 1/v^), = (1/72, Or,-1/72). (23) 

The hadron momenta pi and p 2 have the plus and minus dominant light-cone components, respectively. Accord¬ 
ingly, the quark and gluon momenta ki and i lie along the plus dominant direction while the gluon momentum 
k 2 ~ along the minus direction. The final on-shell photon and quark(anti-quark) momenta can be presented as 


q = yB 


k = xb 



ql 


UbV^ 


* I 

n + q± ■ 


kl 

.n + k± . 


Xb 


7^ 


(24) 
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FIG. 2: The Feynman diagram describing the hadron tensor of the direct photon production. 


The Mandelstam variables for the process and subprocess are defined as 
S=ipi+P2f, T={pi-qf, U={q-p2f, 

s = {xipi + yp 2 f = xiyS, i = {xipi - qf = xiT, u = (q - yp 2 f = yU. (25) 

The amplitude of process ( |22[ ) involves the contributions from 

• the leading (LO) diagrams: two diagrams with a radiation of the photon before and after (^ 2 *^) 

the quark-gluon vertex with the gluon going to the lower blob, see the right side of Fig. 

• the next-to-leading order (NLO) diagrams: eight diagrams constructed from the LO diagrams by insertion 
of all possible radiations of the additional gluon which together with the quark goes to the upper blob, 
see the left side of Fig. 

We denote the sum of diagrams as 

+ ^ 2 ° + ( 26 ) 

that generates the hadron tensor related to the corresponding asymmetry: 

- W = 5151 ' (2'^) 

i=i j=i 

In this work, we mainly dwell on the discussion of the hadron tensor rather than the asymmetry itself. Dia- 
grammatically (Fig. 0, the hadron tensor can be presented in the form of an interference between the LO and 
NLO diagrams: In Fig. I the upper blob determines the matrix element of the twist-3 quark-gluon 

operator while the lower blob - the matrix element of the twist-2 gluon operator related to the unpolarized 
gluon distribution. 


B. Factorization procedure 


Since the collinear factorization is our main tool, let us outline the main stages of factorization. The factor¬ 
ization procedure contains the following steps: 


• the decomposition of loop integration momenta around the corresponding dominant direction: ki = Xip -f 
(ki ■ p)n -\- kx within the certain light cone basis formed by the vectors p and n (in our case, n* and n); 

• the replacement: (Aki <Aki dxi5{xi — ki -n) that introduces the fractions with the appropriated spectral 

properties; 

• the decomposition of the corresponding propagator products around the dominant direction: 


H{k) = H{xp) + 


dH{k) 

dkp 


kl+... 


k—xp 
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• the use of the collinear Ward identity, if it requests by the needed approximation: 


dH{k) 

dkp 


Hp{k, k ); 


• performing of the Fierz decomposition for ippiO) in the corresponding space up to the needed pro¬ 

jections. 

Notice that, for our purposes, it is enough to be limited by the first order of decomposition in the third item. 
As a result of this procedure, we should reach the factorized form for the considered subject, see (©■ 


C. Hadron tensor: QED gauge invariance 


At the first item, we want to discuss the QED gauge invariance of the hadron tensor. To check the QED 
gauge invariance, it is sufficient to consider the typical Feynman diagrams HI, H3 and H5 represented in Fig. 

For the definiteness, we pay our attention on the anti-quark contribution. All our results can trivially be 
extended to the quark contribution as well. 

Before factorization, the Hl-diagram in Fig|^ leads to the following expression: 

W(diag.Hl) = J(d'^ki) 5^'^\ki + fca - fc - g) C 2 x 

v{k)eS{k + J{d*i) S{i + k + q)j1S{k + q)qv{k) ^^2 (28) 

where C 2 implies the corresponding colour factor. Here and in what follows the coupling constants are not 
shown explicitly. In Eqn. ( |28[ ) the gluon (unpolarized) twist-2 parameterizing function is defined as 

^f{k2) I (A2)e-*'=^''= (p2|A“(0)|Px) {Px\A^im)\P2) = gf T%k2) (29) 

while the quark-gluon twist-3 parameterizing function is given by 

I (dV)(d4z)e-*'=i’'^-*^nr,5(pi,ST|V:(r?i)|Px)(Px|^(0)7+A(((z)|ST,Pi). (30) 


We now carry out the standard factorization procedure and, after some algebra, obtain the following expression: 

>V(diag.Hl) =-2 J dxi dy SA\xiPi + yp 2 - k - q) {y) C 2 x 

-7 i ciX 2 - -^— y—71?;— . gu(,/cj«P^ (cci, 0 : 2 ), (31) 


v{k)i 


7 


2xipi +ie J 2 x2Pi + *e 2yp2 + le 


where 

^^l*^’'’{xuX2) = ^ /(rfAi)(dA2)e-“^^^-*("^-"i)^nr,,(pi,ST|V'(Ain)|Px)(Px|^(0)7+A(((A2n)|ST,Pi) 

= e^+‘^^“(pip2)H^(a:i,a:2). (32) 


In the diagrams HI, H3 and H5, we deal with the final state interaction (FSI) with respect to the hard part. 
Therefore, the P^-function has the representation as in (A.131 (cf. [15]), i.e. 


B'2{xi,X2) = B2si{xi,X2) 


T{xi,X2) 
Xi — X2— 


(33) 


In the similar way as in the preceding section, we have to restore the path in the Wilson line with A+-fields 
in the quark-gluon correlators which appear in the diagrams HI, H3 and H5 in Figj^ After straight-forward 








z ] that suggests us to use the representation (A. 131 


calculations, we derive the Wilson line in the form [+oo' 
for our B^-function. 

The contribution of H5-diagram in Fig is equal to zero due to the fact that the photon momentum has the 
dominant minus light-cone component and, therefore, the 7 -structure gives ( 7 ”)^ = 0 . 

We now calculate the hadron tensor term associated with the H3-diagram in Fig|^ it reads 

yy(diag.H3) =-2 J dxi dy S^'^^XiPi + yp 2 - k - q) T^{y)C 2 x 


v{k)i 


7 


2 a;ipi -|- ie 


7 q dx 2 


T 


7 


- 2 yp 2 + ie 2 x 2 Pi + ie 


7Xfc)$y '’^{Xi,X2), 


(34) 


where the B'^-function is also given by the representation (331 or (A.13I. 


After some 7 -algebra, we can check that the contribution of (341 is equal to the contribution of (311 but with 


an opposite sign. Therefore, the sum of all contribution gives us zero (we remind that while the second diagram 
contribution is equal to zero itself) ie. 


>V(diag.Hl) -k W(diag.H3) -f >V(diag.H5) = 0. 


(35) 


In fact, the identity (351 reflects the QED gauge invariance for the hadron tensor. We emphasize that the QED 
gauge invariance takes place owing to the same complex prescriptions in the definitions of B'^-function (Eqn. 
331. In turn, the same complex prescriptions emanate from the final state interaction presented in HI, H3 and 
H5 of Fig|^ Finally, using the similar logical chain as for the DY-process, we can write that 


FSI ^ ^ gauge [-koo , z ] = 1 ^ T{xi,X 2 )^ 

+ le Xi — X2 — le 


QED GI. 


(36) 


We would like to emphasize that the concrete sign of the gluonic pole prescription is not so crucial for the 
QED gauge invariance because here we deal with only one type of interaction which is the final state interaction 
(this distinguishes the case of QCD gauge invariance which is considered below). It is more important to have 
the same prescriptions in all gluonic poles. In other words, from the point of view of contour gauge, we may use 
the same wrong -Vie prescription in the representation of B'^-function for diagrams HI, H3 and H5 depicted in 
Fig|^ But this wrong prescription still leads to the QED gauge invariance. However, after calculation of the 
imaginary parts for the corresponding asymmetry, the wrong prescription definitely plays a negative role. 


D. Hadron tensor: QCD gauge invariance 

We are now in a position to dwell on the QCD gauge invariance of the hadron tensor for the direct photon 
production. To check this invariance, we have to consider four typical diagrams HI, H5, D1 and H9, depicted in 
Figj^ which come from the corresponding ^-process (see, m)- Notice that the gluon enters in the quark-gluon 
correlator as an internal field. For the QCD gauge invariance, we have to assume that all charged particles are 
on-shell, i.e we deal with the physical gluons only. The substantial differences between this case and a pure 
perturbative Compton scattering case are discussed in Appendix B. 

To write down the Ward identity, we need to replace the gluon transverse polarization by the gluon 
longitudinal momentum in the quark-gluon correlator: 

(pi,5^|^(0)7+V'(m)a+(^)l^^,Pi) 

- J£i {pi, S'^\ip{0)'y+-il}{r]i) a+ {£)\S'^,pi), (37) 

where a^{£) stands for the gluon creation operator. The summation over the intermediate states is not shown 
explicitly. Notice that the parametrization of this correlator through B^-function stays unchanged. 

Consider now the contribution of the Hl-diagram in Figj^to the hadron tensor. Before going further, it is 
instructive to begin with the gluon loop integration corresponding to the mentioned diagram. We have 

J{dH)S{£ + k + q)£L{...a+{£)...), 


( 38 ) 
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where we do not explicitly write the operators which are irrelevant at the moment (cf. @). After factorization, 
we obtain 


dX2 


{d'^i)5{x2 — xi — £n)S{£ + k + q)£[,{...a'^ {£)...) = 
j dx 2 S{x 2 Pi + yp 2 ) {x 2 - Xi)pi J {d'^£)S{x 2 - xi- £n){...a~^{£)...) , 


(39) 


where we decompose the hard part around the dominant direction and put £l = (x 2 — xi)pi which is actually 
dictated by the 7 -structure and the momentum conservation. Using all these, we get the following expression: 


yV’(diag.Hl) = — 2 J dxi dy 5^'^\xipi + yp 2 — k — q) {y) C 2 x 

{x2 - a;i)7+7“ 


v{k)e 


1 

2xipi + ie 


7 


dX2 


7 


2 x 2 + 2xipi -|- ie 


—i*v{k) X 


where 


{(-) J{dXi)e (pi,S'^|'0(O)7+'0(Ain) J{d‘^£) 6 {x 2 - xi - £n) {£)\S'^,pi)Y 


(-) {dXi)e (pi, S'-'|V’(O) 7 +' 0 (Ain) {d^£)S{x 2 - xi - £n)a'^{£)\S,pi) = {xi,X 2 ) ■ 


(40) 


(41) 


As can be seen, however, that this diagram does not contribute to the Ward identity. Indeed, after calculation 
of the imaginary part we get the factor (x 2 — xi) in the numerator of (401 which goes to zero owing to d{x 2 —xi) 
from 9mi3''^(xi, X 2 ). 

Further, calculation of the H5-diagram, presented in Figj^ gives us 


^^(diag.HS) = J dxi dy {xipi + yp 2 - k - q)F^{y) C 2 x 


v{k)e ; 


— 

:7 £ 


'dx- 


- v{k) {xi,X 2 ), 


2xipi +ie J 2 x2Pi + ie 
while the contribution of the Dl-diagram in Figj^ takes the form 

>V(diag.Dl) = - j dxi dy S^'^^xipi + yp 2 - k - q)£F^{y) Ci x 

v{.k)e -— . 7 ~ 7 + 7 ~ -— ^ . i* v{k) [ dx2 B^{xi,X2) ■ 

Ixipi + le 2xipi +ie J 


(42) 


2xipi + ie 2xipi + ie 

Finally, the contribution of the H9-diagram with the three-gluon vertex, see Fig|^ reads 

>V(diag.H9) = —4i J dxi dy 6^'^\xipi + yp 2 — k — q) B^{y)C 3 x 


(43) 


v{k)e 


7 


—7 


7 


2xipi + ie 2xipi + ie 


— e* vik) [ dx 2 

J 


X 2 - Xi 


2 (x 2 — xi) + ie 


-B^(xi,X 2 ) . 


(44) 


We now turn to the contour gauge. First of all, based on our previous discussions notice that even a fleeting 
glance is sufficient to anticipate the corresponding prescriptions for i?'-^-functions in (42 H(44l. The H5-diagram 
in Figlifl corresponds to the final state interaction and, therefore, the function BY. should appear here. On the 
other hand, the Dl- and H9-diagrams in Figj^ correspond to the initial state interaction which leads to the 
function B'£^. Performing the explicit calculations (see also m), we can arrive at the conclusion mentioned 
above by restoring the Wilson lines in the quark-gluon correlators of the mentioned diagrams. The Wilson 
line, [-|-oo“, z~], enters in the hadron tensor represented by the H5-diagram in Figwhile the Wilson line, 
[z~, —oo“], appears in the hadron tensor represented by the Dl- and H9-diagrams in Fig|^ 

We sum all contributions and get the following final expression: 


E 


C 2 


N = H1,H5,D1,H9 


W(diag. N) = - —y+y 7+7 y"* 


dx‘ 


BY(Xi,X2) 


X2 


Cl + - + - + 

—27 7 7 7 7 


Sxf 


dX2BY_ (xi,X 2 ) -I- 


4x 


3 + - + 
27 7 7 


' dx - Xi)BY_{xi,X 2 ) 
X 2 — xi+ ie 


( 45 ) 
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where the function BY. are represented by (A.13) or (33), and the function is given by (A.12) or 

BI{xi,X2) = B^Sj^{xi,X2) = . (46) 

xi — X2 + le 

We are now calculating the imaginary part and, ultimately, we derive the QCD Ward identity in the form 

C 2 - Cl - iCg = - , t^] t^e = 0. (47) 


We want to stress that the identity (47) takes place provided only the presence of the different complex pre¬ 
scriptions in gluonic poles dictated by the final or initial state interactions: 

1 r _ . T(a:i,a;2) 


FSI 


ISI 


£+H 

1 


■ le 


-£+ + ie 


gauge [+00 , 2 : ] = 1 ^ 
^ gauge [z~, —oo“] = 1 ^ 


Xi — X 2 — le 
T{xi,X2) 
xi — X2+ ie 


QCD GI. 


(48) 


We emphasize on the principle differences (see details in Appendix 0 between the considered case and the 
proof of the QCD gauge invariance for the perturbative Compton scattering amplitude with the physical gluons 
in the initial and final states. The latter does not need any external condition such as the presence of gluon 
poles. 

Thus, the situation which we discuss is again absolutely similar to that one which has been described in jHlI for 
the dijet production. From (48), it is seen that the different diagrams correspond to the different contour gauges 
and, consequently, to the different functions, B^, which parametrize the hadronic matrix element forming the 
soft part. In this context, we also have a soft breaking of factorization because, first, it spoils the universality 
principle and, second, the gluonic pole prescriptions in the soft part are traced to the causal prescriptions in 
the hard part. At the same time, we can use the replacement |39| : 


T{xi,X2) 
xi — X2 — ie 


T{xi,X2) 
xi — X2 + ie 


27 r iS(xi 


X2)T{xi,X2) 


(49) 


and finally get the same function B^ for all diagrams (in other words, we can use the same contour gauge for 
all diagrams). However, it contains the additional 5{xi — a: 2 )-term which may lead to the collinear factorization 
violation in the same manner as in |39| . The full analysis of this case will be implemented in our forthcoming 
work. 


IV. HADRON TENSOR OF THE DIRECT PHOTON PRODUCTION H: NEW CONTRIBUTIONS 


In this section we calculate the full expression for the hadron tensor which involves both the standard and 
new contributions to the gluon pole terms. The full expression for the hadron tensor related to the case we are 
discussing can be split into two groups: (i) the first type of contributions corresponds to the diagrams H1-H12 
depicted in Fig. [^and, before factorization, takes the following form 


>V(diag.H) = J ( 2 ^) 3 ^^ (2^)32g J{d^ki)id*k 2 )S^^Hki + k 2 - q - k) {k 2 ) x 

y'(d4£) ^(A:i, £) , ^ 2 , £), 

and (ii) the second type of contributions is given by the diagrams D1-D4 in Fig. can be presented as 

W(diag.D) = J (2^)3^^ { 2 tf 2 e J{d‘^ki){d*k2)S^*\ki + k2 - q - k) x 
^fik2)tTD [$«(fcl)^“^(^l,fc2)] . 


(50) 


(51) 


In eqns. (50) and (51), the corresponding coefficient functions are denoted by Hai 3 ^p{ki,k 2 ,£) and Dap{ki,k 2 )- 

v + 1 


The unpolarized twist-2 gluon distribution ^^{k 2 ) and the twist-3 quark distribution 4)^ ^’^{ki,£) are defined 
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presented in Fig. is given by 


in the standard forms, see (29) and (30). The twist-3 quark distribution which appears in the diagrams D1-D4 


“T 'IC J 


(52) 


where the sum over the corresponding intermediate states is implied. We now perform the factorization proce¬ 
dure for Eqns. (50) and (51), and we obtain 


d>V(diag.H) = J +qj-)Cg J dxidySixi - xb) S{y - ys) g'f x 


d^k 


(27r)32F; J (27r)32e 

J dX 2 ^B ^'^{Xi,X 2 ) H°‘>^'P{xi,X 2 ) , 

for the first type of contributions and 

d^q f d^k 


(53) 


dW(diag.D) = 


(27r)32£; J (27r)32e 


d^^\kB+ci±)CD X 


dxidyS{xi - xb) S{y - yB) -g^^{y) 9°^ tr^ [$^^^(a;i) L>“'^(a:i)] 


(54) 


for the second type of contributions. 

To simplify our calculations without losing generality, we may impose the frame where q\ ^ S. The Mandel¬ 
stam variable defined for the subprocess, u, is a small variable and can be neglected. It means that the Bjorken 
fraction yB becomes independent of Xb, and one can write yB = —T/S (due to s -I- t -I- u = 0). 

The next stage is to determine the type of the twist-3 function {xi,X 2 ) which is related to the certain 
complex prescriptions in gluon poles according to the way described in the preceding sections. The correct 
definition of {xi,X 2 )-ixinction should be implemented for each of diagrams. Also, it is instructive to notice 
that the diagrams H1-H8 would not possess the gluon poles in the case where the function B^ assumed to be 
a real one. However, this is not true in our case. 

After computing the corresponding traces and performing simple algebra within the frame we are choosing, 
it turns out that the only nonzero contributions to the hadron tensor come from the diagrams HI, H7, D4 and 
HIO: 


dW(diag.Hl) = 


d^q 


j (2^)32£ ^^^^^^~^ + qj-) C 2 J dxidyS{xi-XB)S{y-yB)J^^{y) x 


(27r)32A J (27r)32e 

2S^ Xi y"^ 




dX2 


[x 2 yS + it][xiyS + ieY vt 


B^{xi,X2) 


(55) 


dW(diag.H7) = 


d^q 


/ (2^)32£ ^^^^^^-^ + q-L)Ci J dxidyS{xi-XB)S{y-yB)J^^{y) x 


(27r)32A J (27r)32e 


[X 2 T + ie] [xiT + ieP 


pi 


(56) 


and 


d>V(diag.D4) = 


{2TTY2E J (2^)32£ ’^^^^*-^-^ / '^3;irfj/i5(a:i-a:B)<5(y-yB)|^®(y) x 

2S'^ xi(y-2yB) 


[xiT + ie 
d^q 


^ 2xiPi + 

d^k 


- dX2BX{xi,X2) , 

le J 


dW(diag.HlO) = J (2^)326 - ^b) Siy - pb) J^’^iy) 


{2Trf2E J (27r)32e 

2T{xi - X2){2T + Sy) 


dX2 


[xiT + ie][x 2 T + ie][[xi - X 2 )yS + ie] 


-g_L+S_L — 

^ dV 


Bi{xi,X2) ■ 


(57) 


(58) 
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Here, Ci = CpNc, C 2 = —Cpl2, C 3 = Cp NcCa/2. The other diagram contributions disappear owing to the 
following reasons: (i) the 7 -algebra gives ( 7 “)^ = 0; (ii) the common pre-factor T + yS goes to zero, (iii) the 
diagrams H2 and H5 cancel each other. 

Analysing the results for the diagrams HI, H7, D4 and HIO (see Eqns. ([5^-( 58 1 ), we can see that 


d>V(diag.Hl) -H dW(diag.H7) -f dW(diag.D4) = dW(diag.HlO). 


(59) 


In other words, as similar to the Drell-Yan process, the new (“non-standard") contributions generated by the 
diagrams HI, H7 and D4 result again in the factor of 2 compared to the “standard" diagram HIO contribution 
to the corresponding hadron tensor. This is our principle result. 


V. CONCLUSIONS 

In this work, we explore both the QED and QCD gauge invariance of the hadron tensor for the direct photon 
production in two hadron collision where one of hadrons is transversely polarized. We present the important 
details related to the use of the contour gauge for the Drell-Yan process with essential transverse polarizations. 
We study the effects which lead to the soft breaking of factorization through the QED and QCD gauge invariance. 

We show that the contour gauges for gluon fields play the crucial role for our study. The contour gauge 
belongs to the class of non-local gauges that depends on the path connecting two points in the correlators. 
It turns out, in the cases which we consider, the prescriptions for the gluonic poles in the twist 3 correlators 
are dictated by the prescriptions in the corresponding hard parts in the similar manner as for the DY-process 
considered in m- 

For the direct photon production, we demonstrate that the prescriptions in the gluonic pole contributions 
differ from each other depending upon the initial or final state interactions in the related diagrams. We stress 
that the different prescriptions are needed to ensure the QCD gauge invariance. This situation has been treated 
as a soft breaking of the universality condition resulting in factorization breaking. Besides, the presence of the 
complex prescriptions in the gluonic pole contributions allow the extra diagrams to contribute nontrivially to 
the hadron tensor. 

We find that the “non-standard" new terms, which exist in the case of the complex twist-3 H'^-function 
with the corresponding prescriptions, do contribute to the hadron tensor exactly as the “standard" term known 
previously. This is another important result of our work. We also observe that this is exactly similar to the 
case of Drell-Yan process studied in na¬ 
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Appendix 


Appendix A: H'^(a;i,*2)—function within the contour gauge 

In this Appendix we give some details regarding the use of the contour gauge. The solution of the QED 
gauge invariance problem for the DY hadron tensor can be found by using the contour gauge conception (see, 
for example, na 1 ^). Here, we would like to demonstrate that there is no an ambiguity in the integral 
representation of the gluon fields through the strength tensor (Eqns. ^ and @). It is important to note that 
the axial gauge condition, A+ = 0, (as well as the Fock-Schwinger gauges) is actually a particular case of the 
most general contour gauge where the Wilson line with an arbitrary path determines the gauge transformations. 
The contour gauge was the subject of very intense studies many years ago. The preponderances of the use of 
the contour gauge is that the quantum gauge theory becomes free from the Gribov ambiguities. 

Let us briefly discuss the main items of the contour gauge conception. To describe this class of gauges, it 
is instructive to assume the geometrical interpretation of gluons where the gluon field is a connection of the 
principal fiber bundle 7 ^(K"^,G, 7 r) (here, implies the base where the principal fiber bundle is determined, G 
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denotes the group defined on the given fiber and tt is a transformation of the base into the fiber bundle V). 
Each element g(a;) of the fiber, with the help of the gluon field Aa, defines the gauge-transformed field: 

= ^~^{x)A^{x)^{x) -f ^g”^(x)a^g(x). (A.l) 

The set of these fields for all g(x) forms the orbit of the gauge-equivalent fields. It is well-known that in order 
to quantize the system of the gauge fields, one has to choose the only element of each orbit. In contrast to 
the usual way, we first fix an arbitrary point (a:o,g(a;o)) ^ in the fiber. Then, we define two directions: one of 
them in the base, the other in the fiber. The direction in the base is nothing else than the tangent vector 
of a curve which goes through the given point Xq- At the same time, the direction in the fiber can be uniquely 
determined as the tangent subspace which is related to the parallel transition. After following this procedure, 
one can uniquely define the point in the fiber bundle. 

Further, solving the parallel transport equation which is defined on the fiber as 

^^VMv))=0, (A.2) 

one can find the solution in terms of the Wilson line: 

g(x) = Pexp^ig J dw • A(w)|g(a;o), (A.3) 

F(xo,x) 


where the points Xq and x are connected by the path P. The starting point Xq is usually fixed, i.e. Xq is 
independent on x. Here, g(a;o) is chosen to be equal to unity Note that the fixing of g(a;) ensures a unique 
choice of the element in the orbit. Inserting (A.3) into (A.l), one can see that the field AJ^{x) is completely 
determined by the form of the path which connects the starting and final points. Moreover, using (A.l) and 
(A.31, one obtains the property: 


X X 

Pexp^ig J du ■ A^{ijj)'^ = g~^{x)Pexp^ig J dw • A(tLi)|g(a;o)- 


(A.4) 


Inserting Eqn. (A.31 into Eqn.( AT| (we remind that the point xq is fixed), we arrive at 

= J dza^g~^{z)Gap{z\A)g{z) = J dza^^ Gap{z\A^), 


(A.5) 


P(xo,a:) 


P(a;o,a:) 


where 


Gcii{z\A^) = Glp{z) =g i(z)Ga/3(2;|A)g(z) 

The contour gauge condition demands that g(a:) is equal to unity for all x belonging to the base, i.e. 


X 

[x, xo] Pexp^ig J dw ■ A(a;)| = 1, Vx 


(A.6) 


(A.7) 


Therefore, within the contour gauge, the field A® (see, (A.5l) becomes 


A-s(x)= J dz^^G^0{z\A), 


(A.8) 


IP(a;o,a^) 


^ We assume that the subspace in the surroundings of an arbitrary point belonging to the fiber can be trivialized. That means we 
can introduce the co-ordinate of the point (x, g(a^)). 

^ The discussion regarding the choice of xq can be found in m- 
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i.e. the gluon field A® is a linear functional of the tensor G^u. Let us briefly comment a choice of the boundary 
conditions for gluons. As is pointed out in jMI, the starting point xq of the path P(a::o,a::) (see, ( A.8[ )) is not 
well defined by construction. In turn, this leads to the presence of th e so- called residual gauge transformation. 
To avoid this ambiguity, it is natural to fix the uncertainty in Eqn. (A.8) by A^(a;o) = 0. In other words, the 
starting point is fixed, and it is independent from the destination of the path P(a;o,a;). 

It is easy to see that the representations (|^ and ( |Io| ) (or the representations ( 12 ) and ( [T^ ) can be derived 
from (A .81 by fixing of the path P(a:o, x) as a straight line connecting the point x with Too, respectively. In 
fact, it means that the representations ( |l^ and correspond to two different contour gauges and give us two 
different representations for two different functions BY{xi,X 2 ) and B ^. {x i, X 2 )- These two representations are 
associated with the final and initial state interactions (Eqns.(|3^ and (|46|)). 

Moreover, in order to get a concrete representation for gluons within the axial gauge, we can explicitly 
parametrize the straight line between the points x and ±oo along the “minus" light-cone direction n~: 


Xais) 


±oo 


= Xa±ric 


lim 

e-fO 


+00 
: 0 


Then, using Eqn. (A.8), one gets 


OO 

(.x) =Tna J ds Ga^iix ± ns) > 


£->■0 


(A.9) 


(A.IO) 


Based on the contour gauge, the representations (A.IOI determine two different contour gauges: A“^(x) cor¬ 
responds to the straight line between x and -|-oo, and A'^~ (x) corresponds to the line between —oo and x. 
Moreover, their projections on the light-cone vector n~ are the same ones: 


n-A““±(a;) = A+ =0. 


(A.ll) 


Thus, we can conclude that the r.h.s. of (12) and (13) correspond to the different functions B^{xi,X 2 )' 
BX{xi,X 2 ) =S{xi- X2 )^I(_oo)(2;i) + for the gauge 

Xi X2 \ '26 

T{xi,X2) 


BI{xi,X2) = 5 {xi - X2)B\f^^^~^{xi) + —- 


xi — X2 — le 


for the gauge 


[a:, -oo] = 1, 
[+00, x] = 1 . 


(A.12) 
(A. 13) 


It is clear that these two representations can not be equivalent of each other. Also, it is reasonable to assume 
the zeroth boundary conditions for gluons to be = 0 which is in agreement with Refs. [151136| . Note 

that the functions B^ do not possess a certain property under the time-reversal transformation. This property 
becomes a well-defined one only after calculation of the imaginary part for the hadron tensor. 

We thus have no the ambiguity in the solutions of (|^. Asa result, the R^-function has a non-trivial imaginary 
part which contributes to needed for the gauge-invariant set. 


Appendix B: Comparison with the perturbative Compton scatering amplitude 

In this Appendix, we discuss an essential difference between the QCD gauge invariance for the “perturba¬ 
tive" Comptom amplitude, {q{k)g{£)\S\g{k 2 )q{ki)), and the “nonperturbative" analog of Compton amplitude, 
{X{Px)q{k) \ S \g{k 2 )A{P)), where one of gluons is included in the loop integration (see, Figsj4][§. 


1. QCD gauge invariance of the perturbative Compton amplitude 

We first consider the “perturbative" Compton amplitude represented in Fig|^ To proof the QCD gauge 
invariance, we assume that all gluons are physical ones with the transverse polarizations. The first diagram in 
Fig|^ gives us 


AapidiaA) = {-i)gu{k)^p 


ki + k2 


{ki + ^ 2)2 -I- ie 


7 q, P t u{ki ), A(dia.l) = Aap{dia.A)e-^ Pp . 


(B.l) 
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To check the gauge invariance with respect to, for example, the physical gluon with momentum I, we h ave t o 
replace the given polarization vector on the momentum provided I ■ (* ^ = 0, in the amplitude (B.l). 


We have 


yl(dia.l) = {—i)gu{k)i 


ki + k2 


(fci + ^ 2 )^ + ie 


u{ki). 


(B.2) 


Using the momentum conservation, ki + k 2 = k + £, and the corresponding equations of motion, the Eqn. (B.21 
reduces to the following form: 


^(dia.l) = {—i)gu{k) [fci +k 2 +k] 


‘ u{ki) = uik) e"*' u{ki). 


{ki + ^ 2)2 + it 

In the similar way, we can get the contribution of the second diagram presented in Fig|^ It reads 

ki —. 


^(dia.2) = {—i)gu{k) e'^ 


-I t” u(fci) = — u{k) t'^ u(fci). 


(fci — £Y + it 

Therefore, the sum of these two diagrams gives us the color commutator: 

^(dia.l) + ^(dia.2) = {—i)gu{k) t'^ [t“, t^] u{ki). 


(B.3) 


(B.4) 


(B.5) 


In the Feynman gauge^ , the third diagram in Fig. 0 contributes as (where we again replace the transverse 
gluon polarization on the gluon momentum: i) 


^(dia.3) = igu{k)"ffjt‘^ u{ki) 


(k — fci)^ + it 


~ir^V^ap{-k2;k-ku£)tt£0.. 


(B.6) 


where 


kaa'/3( ^2, k ki^ — gatjik ki T fe2)/9 k ki')Q> gal3{_k2 “t“ £)a • (^■^) 

Making use of e-*- • ^ = e* = 0 and 


{k — ki + k2) ■ £ _ k2 ■ £ _ 

(k-ki)^ “ (fe - £)2 " " ’ 

we obtain 

M(dia.3) = i^g f°‘^‘^u{k)t'^ u{ki). 
The sum of all three diagrams leads to 

yl(dia.l) + ^(dia.2) + yl(dia.3) = 0 


(B.8) 


(B.9) 


(B.IO) 


which ensures the QCD gauge invariance. In the case of the “perturbative" Compton amplitude with physical 
gluons, the momentum conservation, S^‘^\ki + k 2 — £ — k), and the on-shellness of all external particles play 
important role for checking of the gauge invariance. 


2. QCD gauge invariance of the nonperturbative Compton scattering amplitude 

In order to get the diagrams in Fig. we now attach the nonperturbative blobs to the diagrams in Fig. 

In contrast to the “perturbative" Compton amplitude, the gluon with momentum, £, with respect to which 
we check the gauge invariance, is included in the loop integration. As a result, the momentum conservation for 
the subprocess has a form: ki + k 2 = k. 


® For this discussion, the type of gauge is really not important. 
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Let us focus on two diagrams in Fig. which are enough to demonstrate how the color commutator contri¬ 
bution can be formed. Both diagrams are given by 


{X{Px)qik)\S\g{k 2 )A{P)) = u(fc)e„ / (d ^){d r])e 




{Px\ 




S'4){ri)dAai^) 


\P) 


Ij>=...-A=0 


The first diagram in Fig. [^contributes as 

S(dia.l) = u(fc) ({d^tj^pS{i + k)e^ ({d^vi)e~''^'^{Px\Ap{vi)t°'t'"'ij;{Q)\P) . 


(B.ll) 

(B.12) 


This expression corresponds to the amplitude before factorization and, therefore, the quark with ki + ^ and 
gluon with ^ are off-shell. Next, we apply the factorization procedure and derive the following expression for 
the first diagram contribution (here we use the light-cone basis presented above): 


e(dia.l) = (B.13) 

u[k) j dx 2 {x 2 — X\)n*S{^{x 2 — xi)n* + k)i^ P u{x 2 n*){Px\a'^ {{x 2 — Xi)n*')b~ {x2n*)\P) 


or 


S(dia.l) = u(A:) J dx2"f P u{x2n*){Px\aA {^{x2 — Xi)n*)b {x2n*)\P) , (B.14) 


where we replace the transverse gluon polarization vector by the gluon momentum (x2 — xi)n* . In Eqn. (B.14), 
we use the Fourier transformations for the quark and gluon fields, a"*" and b~ denote the gluon creation and 
quark annihilation operators, respectively. 

In the similar manner, we consider the second diagram contribution of Fig. [^ Before factorization, we have 


B{dm.2) = u{k)e^S{ki) / P u{ki + ^){Px\a+{^)b-{ki + t)\P). 


(B.15) 


Having applied the factorization procedure, we derive the following expression: 


H(dia.2) = ^(A:) J dx 2 e'^j P u{x 2 n*){Px\a~'~ [{x 2 — xi)n*)b {xin*)\P) . (B.16) 


Analyzing Eqns. ( |B.14 | and (B.161, one can see that in order to form the color commutator combination we 
have to insist on the “external" conditions which actually emanate from the presence of gluon poles. Indeed, if 
the gluon poles are present, the amplitude H(dia.l) is accompanied by the factor of iTrS{x 2 — xi). At the same 
time, the amplitude H(dia.2) has the factor of —iTT 6 {x 2 — Xi) in according to the FSI and ISI prescriptions, see 

(|4§. 

Thus, to check the QCD gauge invariance, the case of the “perturbative" Compton amplitude with the physical 
gluons in the initial and final states does not need the external condition for the presence of gluon poles. It is 
sufficient to use only the momentum conservation, ki + k 2 = i + k, and the equations of motion for the initial 
and final quarks. At the same time, in order to demonstrate the QCD gauge invariance for the case represented 
in Fig. where one of gluon momenta belongs to the loop integration, the existence of the gluon poles with the 
corresponding FSI and ISI prescriptions must be included. 
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FIG. 3: The typical Feynman diagrams, contributing to the hadron tensor. 


























































































































































































